Electronic structure and magnetism for FeSi(i_ x )Ge a; from supercell calculations. 
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Recent studies of FeSiji-^jGea;, which found a transition 
from an insulating to a magnetic metallic state near x=0.25, 
have revived the discussion about the role of strong correla- 
tion in these systems. Here are spin polarized band calcu- 
lations made for 64-atom supercells of FeSifi-^jGe^ for dif- 
ferent x and different volumes for large x. The results show 
that the small band gap in FeSi is closed for x > 0.3, be- 
cause of both substitutional disorder and increased volume. 
Ferromagnetism appears near this composition and becomes 
enforced for increasing x. The ^-dependence of the electronic 
specific heat can be understood from the exchange splitting 
of the density-of-states near the gap. Strong volume depen- 
dencies for the properties of FeGe suggest experiments using 
pressure instead of x for investigations of the gap. 
PACS: 75.10.Lp, 71.23-k, 71.30+h, 75.50.Bp. 



I. INTRODUCTION 



The unusual properties of the cubic compound e-FeSi, 
investigated early by Jaccarino et al 1 , continue to be the 
subject of many theoretical and experimental studies 2 - 15 . 
Band calculations based on the local density approxima- 
tion (LDA) 16 , showing that FeSi is a very narrow band 
gap semiconductor with a gap of about 6 mRy 9 - 15 , can 
explain many observed low-temperature properties. The 
sharpness of the density-of-states (DOS) below the gap is 
observed as a dispersionless feature in the spectra of an- 
gular resolved photoemission 4 taken at low temperature, 
T. Optical measurements reveal that the gap is gradually 
filled at higher T, and an unusual feature appears to be 
that the spectral weight is not conserved 5,6 . The problem 
is also to understand why the material behaves as a metal 
with large magnetic susceptibility x(T) at higher T. Elec- 
tronic excitations within the band structure (determined 
at T=0), described by the Fermi-Dirac function, are not 
able to explain such T-dependent properties 1,2 . Many 
theories, mostly based on some form of strong correla- 
tion, have been proposed for an explanation of the un- 
usual T-dependence 7,8,10,11 . On the other hand, elec- 
tronic structure calculations within the LDA can explain 
the T-dependent properties when the effect of thermal 
disorder is included 13 . The filling of the gap is an effect 
of structural disorder for T larger than about 150 K, and 
the material behaves as a Stoner enhanced paramagnet 
above this temperature 13 . 

Different dopings of FeSi or other compositions like 
FeGe, MnSi, etc. lead to a large variety of properties, and 
metallic magnetism is often found 17 . The LDA approach 



has been used to analyze thermoelectric properties in 
doped FeSi 18 . Recent studies 19,20 focused on Ge substi- 
tutions on Si sites, in which the evolution of the gap and 
other properties were followed continuously between the 
isoelectronic FeSi and FeGe systems. Pure FeGe is metal- 
lic and magnetic with a long-range spiral spin alignment, 
and the saturation moment is about 1 /j,b per Fe 19,21,22 . 
The alloy FeSi( 1 _ 2 ,)Ge x shows a transition from an insu- 
lating state to a metallic ferromagnetic one near £=0.25, 
with enforced magnetism at larger x 19 . Some, but not 
all properties are interpreted in terms of the Kondo insu- 
lator model with a Hubbard U as parameter 19 . In order 
to pursue the search for the alternative solution based on 
LDA band structures, we study here the evolution of the 
gap and the magnetic moment as function of Ge substitu- 
tions with special attention to the effect of substitutional 
disorder. Methods like the virtual crystal approxima- 
tion, which interpolate potential parameters between the 
pure constituents, are insufficient to study the effect of 
disorder. Therefore we use supercells with different site 
occupations for the band calculations of the FeSi(i_a;)Ge x 
alloys. 



II. METHOD OF CALCULATION. 



The bandstructure for the ordered material can be de- 
scribed by the use of a supercell, although it is a more 
complicated calculation than for the normal cell. The 
bands of the normal cell are folded back into the smaller 
Brillouin Zone of the supercell, giving more bands per 
k-point. The bands are degenerate along different direc- 
tions in k-space, but if there is disorder (structural or 
substitutional) there will be slighly different band dis- 
persions along different directions. The average energy 
of a band along say k x and k y may be similar to the en- 
ergy for the ordered case, but the non-degeneracy along 
different directions can be interpreted as a band broad- 
ening of the original band. Thus, the band broadening 
or smearing of the bands near the gap, which becomes 
important for the properties of the FeSi-like materials, is 
because of a number of non-degenerate states that devi- 
ate from the the original energy as the degree of disorder 
increases. These effects, as well as a possible drift of the 
average energy due to disorder, are contained in the su- 
percell approach. Disorder of real alloys does not show 
the long-range periodicity of the supercell, but this a mi- 
nor problem when large supercells can be used. 
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The spin-polarized version of LDA can be used to de- 
termine the magnetic moment, but the question of spiral 
order of the spin moments in FeGe is not addressed in 
this work. The calculated moments will be compared 
with saturation moments at high field (0.3 T) at which 
the moments are aligned. Since a magnetic field of this 
amplitude is very small on the energy scales in a band 
calculation (1 mRy corresponds to 230 T) it is unlikely 
that this field can modifiy exchange splitting or the elec- 
tronic structure. 

The self-consistent Linear Muffin-Tin Orbital method 
is applied to 64-atom supercells (8 basic unit cells of 
the B20-structure) as is described for the calculation for 
structurally disordered FeSi 13 . The linearization energies 
are taken near the center of each ^-band. The LDA po- 
tential contains no special on-site correction due to corre- 
lation. All sites are nonequivalent in the disordered cells. 
The self-consistent iterations use initially 8 k-points in 
1/8 of the Brillouin Zone, while they are terminated by 
iterations using 27 k-points. The band gaps in the two 
sets of k-points agree to within 0.5 mRy, but the mag- 
netic moments in the metallic cases can differ by 20-30 
percent for some configurations. This difference comes 
from difference in the DOS near the gap between the 
two sets of k-points. A comparison of the DOS from a 
case using of 27 and 64 k-points show no large difference 
and it is expected that the convergence of the moments 
from 27 k-points is satisfactory. The calculations include 
a thermal smearing of 1 mRy in the electronic occupa- 
tion (Fcrmi-Dirac) . This stabilizes the self-consistency, 
but there is also a physical reason behind the smearing; 
The zero-point motion of the atomic positions, leading to 
some structural disorder, can be estimated to be of the 
order 0.5 percent of interatomic distances for appropriate 
values of force constants and atomic masses 2 ' 15 . From 
various arguments it can be concluded that an essential 
effect of band smearing on the electronic structure near 
the Fermi energy, Ep, can to some extent be modeled 
through the Fermi-Dirac occupation. Although the role 
of structural disorder (zero-point motion, thermal disor- 
der as function of temperature and possible disorder be- 
cause of Si/Ge substitutions) merits a more careful anal- 
yse, we may include some of the effects from zero-point 
motion through the Fermi-Dirac function. All calcula- 
tions are made with the same internal structural parame- 
ters of the B20 structure and with no structural disorder. 
Thus, only the effect of substitutional disorder appears in 
these calculations, although it is probable that different 
relaxation around Si and Ge will increase the effects of 
disorder in real FeSi(i_ 2; )Ge a ;. 

The calculations do not determine the equilibrium vol- 
umes through minimization of the total energies for each 
composition, but the calculated pressures (P) are used as 
a guidance for finding the volumes relative to that of pure 
FeSi. Absolute values of direct calculations of the pres- 
sure are less reliable than when P is calculated through 
the volume derivative of total energies, but they converge 
more rapidly and can be useful for studies of relative 



variations. The bands of pure FeSi at T=0 (no disorder) 
have a gap of 6 mRy at the Fermi energy, E^, in agree- 
ment with other calculations using different methods 9 - 
14 . This is for a lattice constant, a =4.39 A, between 
the experimental one (4.52 A) and the theoretical one 
(at the minimum of the total energy). This reflects the 
usual problem of using the LDA for 3d transition metals, 
in which the lattice constants come out 2-3 percent too 
small compared to experiment 23 . This uncertainty could 
be relatively severe for FeSi because of the narrow gap. 
However, the theoretical E 9 agrees well with experiment 
and it decreases only slowly when the lattice constant 
increases, in agreement with the pressure dependence of 
the magnetic susceptibility 12 ' 14 . 

The first calculation for FeGe is made with a =4.52 
A, i.e. about 3 percent larger than for FeSi. This choice 
is based on the differences in covalent radii and lattice 
constants of pure Ge and Si, and as for FeSi we con- 
sider the LDA-bands for a lattice constant which is a 
few percent smaller than the experimental value (4.7 A). 
The calculated gap in FeGe, 1.5 mRy, is close to the 
LDA result of ref 20 . For intermediate compositions a 
is assumed to vary linearly with x, which turns out to 
be consistent with the calculated variations of the pres- 
sure. The discontinuities of the potential at the limits 
of the different Wigner-Seitz spheres are small when the 
Wigner-Seitz radii are 1.44 A for Fe, and 1.37 A for Ge 
and 1.27 A for Si in FeGe and FeSi respectively. Almost 
the same values are maintained in the calculations for the 
FeSi( 1 _ a; )Ge a ;-systems. The s- and p-bands of Si and Ge 
are similar, both concerning the position relative to the 
Fe bands and the band widths. The unoccupied 3d band 
in Si is almost 0.5 Ry closer to Ef than the 4d band in 
Ge. This makes a difference for the hybridization with 
the Fe-d bands near and for the degree of disorder of 
the band structure in the alloys. 

Calculations are made for pure Fe32Si32 and Fe32Ge32 
and for configurations with 4, 7, 10, 11, 12, 14, 16, 17, 
19, 22, 25, 28, 30 and 31 number of Ge sites per super- 
cell. The distribution of Ge vs. Si sites is random in 
the different cells for different x, but two calculations are 
made in which the Ge sites (4 and 12, respectively) are 
clustered together within the cell. Two calculations are 
made for pure FeGe, labeled II and III, in which ao is in- 
creased by one and two percent, respectively. As will be 
dicussed later, the results for the moment agree best with 
experiment in calculation FeGe-III, although the lattice 
constant is smaller than the experimental one. 
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III. RESULTS. 



A. Bandstructure near the gap. 

The width of the gap, the unpolarized DOS at E^, 
Npara, and the exchange splittings, £, from the spinpo- 
larized calculations are shown in Table I. The gap, E g , 
is about 2.5 mRy in both calculations with 4 Ge atoms, 
either randomly distributed or clustered within one basic 
cell of FeGe together with 7 basic cells of FeSi. The effect 
on the gap from Ge substitutions is relatively weak com- 
pared to other dopants replacing Fe or Si 18 . As might be 
expected, the effects of disorder on the DOS are to reduce 
Eg and to make the DOS peaks near the gap wider. The 
gap is reduced when the concentration of Ge is increased, 
and the alloy becomes metallic at about 8-9 Ge sites per 
cell. The Fermi energy is near the bottom of a 'valley', or 
pseudogap, in the DOS for all paramagnetic calculations, 
showing that the gap is never completely washed out. A 
real gap is restored for the two highest Ge concentrations 
when only one or two Si remain in the cell. 

The Fe d-bands contribute most to the large DOS near 
the gap, with some local variations from site to site due to 
the differences of near neighbor atoms. Magnetic, Stoner- 
type ordering can start at sites with the largest N(Ep). 
The magnetic instability occurs when N(Ep) ■ I > 1, 
where / is the exchange integral, but the size of the mo- 
ments is determined by other parameters. (As seen in 
Table I, the calculation with 12 clustered Ge atoms has 
a higher paramagnetic N(Ep) than the case with ran- 
dom site occupation, but for £ and the moment it is the 
other way around.) Typically there is more than a factor 
of two between the lowest and highest local moments on 
different Fe when disorder plays a role, i.e. for x not too 
close to 1. The difference in local moments leads to an 
additional disorder in the spin-polarized part of the po- 
tential, so that the majority and minority DOS functions 
are not exactly like a rigid splitting of the paramagnetic 
DOS. This makes the remaining pseudogaps in the DOS 
of the two spins more 'washed out' than in the param- 
agnetic DOS. The mechanism leading to magnetism is 
accelerated as soon as the exchange splitting allows for 
Ep to enter into the high peaks of the DOS of majority 
and minority spin above and below the gap, respectively. 
This process slows down the self-consistent convergence 
in some cases. As will be shown later, it also provides 
an explanation for the large variation of the moment in 
FeGe as function of volume. 

The spinpolarized N(Ep) and magnetic moments m 
are summarized in Figs. 2 and 3. Calculations for 4 and 
7 Ge per cell give vanishing moments, while those for 
30 and 31 Ge have finite moments despite a small gap. 
The latter can be understood from the thermal smearing 
that overcomes gaps smaller than about 1.5 mRy. These 
results show that weak magnetism can be the result of 
substitutional disorder of the FeSi(i_ :c )Ge a: -system for x 



larger than ~ 0.3. This range for ferromagnetism agrees 
well with experiment 19 , but the moments are smaller, 
and as will be discussed later, the moments are sensitive 
to volume. 



B. Specific heat. 



The DOS near the gap is very peaked, and statisti- 
cal fluctuations of N(Ep) can be seen for the different 
configurations. The total, spin-polarized N(Ep) values 
shown in Fig. 2 are scattered around 800 states/Ry/cell 
for large Ge concentrations, while the maximal values 
(1000-1200 states/Ry/cell) are found for x * 0.5. The 
electron-phonon coupling A is calculated to be about 0.2 
for disordered FeSi having a DOS of 700 states/Ry/cell 13 . 
As A is proportional to N(Ep) one can estimate the elec- 
tronic specific heat coefficient 7 = ^n 2 k B N(Ep)(l + A) 
to be « 5-9 mJ/mole K 2 when N(Ep) varies between 800 
and 1200 states/Ry/cell. This is lower than measured 19 . 
Such discrepancies are often attributed to spin fluctua- 
tions that contribute to 7 through X s f in near magnetic 
systems. A calculation of X s f for a single 8-atom cell 
of disordered FeSi (which becomes metallic with a mod- 
erate DOS below the limit for Stoner magnetism) indi- 
cates that it can be large, between 0.5 and 1. This DOS 
is comparable or lower than the DOS in FM phases of 
FeSi a; Ge( 1 _ a; ), but X s f is often lower on the magnetic side 
of a FM transition 26 . Therefore it can be expected that 
spinfluctuations contribute to 7 just at the very begin- 
ning of FM and less when magnetism becomes stable at 
larger x. The measured 7 is peaked at ~ 20 mJ/mole 
K 2 near x = 0.37 19 , i.e. rather close to the critical x 
for the metallic transition. From the DOS of Fig 1 it is 
possible to understand this behavior. All DOS functions 
show the peaks around the gap or pseudogap indepen- 
dently of x. Larger disorder will smear the DOS, but the 
peak positions remain even in the spin polarized major- 
ity and minority DOS functions. An inspection of the 
DOS functions shows that a moderately large moment 
of 3 [1b per 64 atom cell (corresponding to a £ of about 
9 mRy) will have Ep on the first peak above the gap 
within the majority DOS, and on the second peak below 
the gap within the minority DOS, cf. Fig 1. Thus the 
total N(Ep) is large for a moment of this size, from Fig. 
1 one can estimate that the combined N(Ep) (majority 
plus minority) is about 1200 states/cell/Ry. If the mo- 
ment is larger, then Ep is found further to the right on 
the majority DOS and further to the left on the minor- 
ity DOS. N(Ep) are lower for both spins, and the com- 
bined N(E F ) taken from Fig. 1 is 700-800 states/cell/Ry 
for a large moment (7-8 hb /cell) and a large £ of 12-15 
mRy. This discussion can be exemplified by rigid-band 
spin splittings on the paramagnetic DOS, as in Fig. 4 
for three different compositions. The total charge is con- 
served, whereas a moment W is obtained from the dif- 
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ferences in the number of occupied spins. If m shows a 
linear increase as function of x it is expected that "f(x) 
has the same shape as in Fig. 4. 

The trend of increasing moments towards large Ge con- 
centration is evident from the points in Fig. 3. The cal- 
culated moments are typically 2-3.5 /i B for x in the range 
0.35-0.7, while for larger Ge concentrations the moments 
can reach 5 /(b or more. Thus, the former range of con- 
centration corresponds to a large N(Ep) (and 7), while 
in the latter N(Ep) is reduced. (The relative variations 
in 7 will be enhanced over those in N{E F ) because of 
the coupling factors A.) This behavior is partly con- 
firmed in Fig. 2; N(Ep) is largest for x between 0.4 
and 0.6, although there is scattering among the points 
due to the small set of configurations. The relative vari- 
ations of m(x) and 7(2;) fit reasonably well to the mea- 
surements by Yeo et al, but the moments are too small 
when x — > 1. The measured moment at x=0.37, about 
0.1 jUe/f.u. (equal to <~ 3.2[i B for the 64 atom cell), is in 
fair agreement with the calculations. But the moments 
for x — ► 1, ^0.4 fj, B /f.u. (<~ 12fJ, B per 64 atom cell) 19 , 
are much larger than in the calculations, and they are 
not yet the saturation moments. 



C. Properties at Ge-rich compositions. 

The calculations for x — > 1 indicate that a small gap 
will reappear in FeGe if spin-polarization could be pre- 
vented. Calculations for undoped FeGe are easier as they 
can start from results for the smaller 8-atom cells, and 
a few volumes have been studied. The calculated pres- 
sure for FeGe (FeGe-I at a =4.52 A) is about 0.1 Mbar 
larger than for FeSi (at ao=4.39 A), which gives a hint 
that the calculations for x — > 1 are made at too small lat- 
tice constants. In a calculation (labeled FeGe-II in Table 
I), in which ao is increased to 4.56 A , the difference in 
pressure is reduced to about 0.05 MBar, and for FeGe-III 
when ao is 4.61 A the pressure is almost the same as for 
FeSi. The (paramagnetic) gap becomes smaller as a is 
increased (see Table I) . It is unusual that a gap becomes 
wider with applied pressure. It is more common that the 
gaps are between different sub-bands, so that the gaps 
become narrower when the bands become wider at larger 
pressure. But the gap in FeSi and FeGe is within the Fe- 
d band, so the gap widens together with the band when 
the lattice constant is reduced 12 , although the effect is 
small for FeGe as seen in Table 1. A narrower gap for in- 
creased lattice constant implies that a metallic transition 
is approaching and the moment in the spin-polarized cal- 
culations increases from 0.08 /xs/f.u. in FeGe-I to 0.25 
and 1.06 /j, B /i.u. for FeGe-II and III, respectively. The 
increase in moment as ao is increased appears very large 
when there only is a small reduction of the gap. However, 
in addition to a narrower gap there are sharper increases 
of the DOS on both sides of the gap when the volume is 



increased, so that more states come closer to Ep. This is 
probably more relevant for the evolution of the magnetic 
moment as function of volume than the value of the gap 
itself. 

The strong volume dependence of m is extended 
towards increasing Si-compositions, but not too far. 
When the lattice constant is increased by 1 percent for 
Fe 32 Sii4Gei8 and Fe32Si7Ge 2 5, it is found that the mo- 
ment increases only in the latter case (by more than 50 
percent) compared to the moments shown in Figure 3 at 
the respective composition. Thus, when the gap is just 
closed and the system is not yet a good metal as in the 
former case, there is not a strong dependence on volume. 
The DOS has 'tails' near the gap, partly due to disorder. 
The initially small overlap between the DOS of the con- 
duction band and the valence band makes the system a 
rather poor metal, and magnetism is still hesitant. But 
in the latter case, when £=0.78 (as well as when x=l) 
the DOS is large or increases rapidly near a narrow gap 
and the effect of pressure is large. 

As was discussed above, this dramatic increase of m is 
possible when the exchange splitting is sufficiently large 
to make the system clearly metallic within both spins. It 
might seem strange that a moment develops depite the 
(small) gap in the DOS for Ge-rich compositions, and it 
could be suspected to be a metastable state. However, 
convergence of the total energies for FeGe-II show that 
the magnetic state is indeed the stable one with a total 
energy more than 30 mRy lower than the non-magnetic 
one. The interesting situation of a stable non-magnetic 
state and a metastable magnetic state with larger to- 
tal energy is more probable for larger gaps, and the T- 
dependences of the two states could be very different so 
that metamagnetism, large magnetoresistance and phase 
transitions can be imagined. 

The moment is very sensitive to volume and disor- 
der, but the combined results of the pressure calculations 
(compared to that of FeSi) and the size of the calculated 
moments (compared to the measured one) indicate that 
the lattice constant should be close to the one used for 
the calculation FeGe-III, or near 4.6 A, for an optimal 
description of the properties of FeGe. The calculated 
moment is then close to the saturation moment found 
for fields larger than 0.3 T 19 ' 21 - 22 , of the order 32 \i B per 
supercell. It also follows that magnetism on the Ge-rich 
side of FeSi(i_ a; )Ge a ; depends more on the increased vol- 
ume than on disorder. The strong volume dependence 
of m suggests that thermal expansion, apart from effects 
due to thermal disorder, should lead to an unusual in- 
crease of m with temperature. This hypothesis is cor- 
roborated by the measured increase of m(T) by about 10 
percent between low T and just below the Curie temper- 
ature (Tq ~ 280 K) at which m drops to zero, seen in 
the data for x = 1 by Yeo et al. By using a typical coef- 
ficient for thermal expansion as for Fe, one arrives at an 
increase of the lattice parameter of the order 0.2 percent 
between low-T and room temperature. The calculated 
increase of m as function of a translates into an even 
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larger increase of m(T) within this temperature interval, 
but this does not take into account the type of moment 
disorder that finally leads to zero effective moment at Tq- 



IV. CONCLUSION. 



tal information about the existence of a small gap (1-2 
mRy) in the spin polarized bands of pure FeGe. There- 
fore, experimental results of FeGe as function of pressure 
would be valuable. 
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In conclusion, it has been shown that LDA band 
calculations can give an adequate description of the 
FeSi( 1 _ a; )Ge a ; system if disorder and changes in volume as 
function of x are accounted for. The gap is small enough, 
so that the combined effects of disorder and increased vol- 
ume lead to zero gap for x « 0.3. The system becomes 
magnetic for larger Ge concentrations, and the behav- 
ior of the electronic specific heat is directly related to 
the DOS structures of the majority and minority bands. 
The quantitative agreement with experiment for the mo- 
ments is best when the lattice constants for Ge-rich com- 
positions are 1-2 percent larger than in the results shown 
in Fig. 3. 

Pressure experiments can be suggested from the calcu- 
lated results of pure FeGe at different volume. The del- 
icate balance between magnetism in metallic FeGe and 
absence of magnetism in semi-conducting FeGe could be 
followed continuously as function of increasing pressure. 
The effect of uniform pressures would be the same as 
varying x from FeGe towards FeSi, but with the advan- 
tage that effects of substitutional disorder can be ex- 
cluded. Theoretical estimates of the bulk modulus (B) 
within LDA are often larger than what is found exper- 
imentally, as is the case for FeSi 12 . From the results 
above, leading to B w 1.6 Mbar, it is expected that a 
pressure of 0.1 Mbar will be sufficient for a suppression 
of the moment. This represents an upper limit, since 
the measured low value of the Dcbye temperature for 
FeGe 19 is an indication of a low B. Information about 
disorder, also structural and thermal ones, is important 
for calculations of the properties of the isoelectronic al- 
loys FeSi( 1 _ :E )Ge a ; because of the large DOS peaks close 
to a very small gap. 

Finally, the fact that properties depend on a tiny gap 
between high DOS peaks motivates a comment about 
results of band calculations. When the DOS near the 
Fermi energy is rather flat, as in most metallic materials, 
there are no big consequences of details (choice of basis, 
linearization energy, size of atomic spheres, general po- 
tential, type of density functional and so on) in the band 
theory method. But here for FeSi and FeGe, when the 
DOS varies from zero to very large amplitude within 1-2 
mRy, it might be that differences in the method of calcu- 
lation lead to quite different properties. It is important 
to note that the gap of about 6 mRy in FeSi agree with 
experiments and other calculations. Also the result that 
FeSi( 1 _ a ,)Gc a; becomes metallic for x larger than about 
0.3 agrees with experiment. But there is little experimen- 
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TABLE I. Calculated band gap and density-of-states at 
Ep, N para {Ep) 1 in units of states/cell/Ry from param- 
agnetic calculations, and exchange splitting, £, from the 
spin-polarized calculations. The labels "-c" mean clus- 
tered configuration (see text), and I, II and III for pure 
FeGe are for lattice constants 4.52, 4.56 and 4.61 A. 



Cell 


E 9 (mRy) 


Npara 


£ (ml 


Fe32Sl32 
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Fe 32 Si28Ge 4 


2.5 
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Fe32Si2sGe4-c 
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Fe 3 2Si 2 5Ge 7 


2 
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Fe 32 Si22Geio 


- 


16 
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Fe 3 2Si2iGeii 
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Fe 3 2Si 2 oGei2 
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110 
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Fe32Si2oGei2-c 
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Fe 3 2Sii 5 Gei 7 
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110 
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27 


12 


Fe 3 2Si 4 Ge28 
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13 


Fe 3 2Si 2 Ge 3 o 


0.5 





12 


Fe 3 2SiiGe 3 i 
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18 
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FIG. 1. Paramagnetic density-of-states of FeSi, FeGe and 
FeSio.56Geo.44 as calculated for the supercells containing to- 
tally 64 atoms. The energy is relative to Ef- 
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FIG. 3. Magnetic moments as function of the Ge concentra- 
tion x in cells containing 64 atoms. Note that the two calcu- 
lations for pure FeGe (x=l) for 1 and 2 percent larger lattice 
constants have much larger moments, 8.2 and 33.8 /j,b /cell re- 
spectively. Calculations for 1 percent larger lattice constants 
for £=0.44 and :r=0.78 give increased moment only in the 
latter case, near 8 fin /cell. 
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FIG. 2. Spinpolarized density-of-states at the Fermi en- 
ergy, N(Ef), as function of the Ge concentration x in unit- 
cells containing 64 atoms. 
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FIG. 4. Total N(Ef) as function of magnetic moments cal- 
culated from a rigid spin splitting of the paramagetic DOS 
functions in Fig. 1. The heavy line is obtained when using 
the DOS for x~0, the broken line for £=0.44 and the thin 
line for x=l. The variations of the observed specific heat of 
Yeo et al (shown in the inset of fig. 4 of ref. 19 ) show one peak 
near :r=0.37, corresponding to a small m. 
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